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Abstract. The propagation of waves in soft dielectric elastomer layers is investigated. To this 
end incremental motions superimposed on homogeneous finite deformations induced by bias electric 
fields and pre-stretch are determined. First we examine the case of mechanically traction-free layer, 
which is an extension of the Rayleigh-Lamb problem in the purely elastic case. Two other loading 
configurations are accounted for too. Subsequently, numerical examples for the dispersion relations 
are evaluated for a dielectric solid governed by an augmented neo-Hookean strain energy. It is found 
that the the phase speeds and frequencies strongly depend on the electric excitation and pre-stretch. 
These findings lend themselves at the possibility of controlling the propagation velocity as well as 
filtering particular frequencies with suitable choices of the electric bias field. 

Keywords: dielectric elastomers; electroelastic waves; finite deformations; soft actuators; non-linear 
electroelasticity. 

1. Introduction 

The goal of this work is to investigate the propagation of electromechanical induced waves in a 
dielectric elastomer (DE) layer subjected to finite deformations. In the purely elastic case the first 
solution for surface waves based on the exact equations of (2D) elasticity was introduced by Lord 
Rayleigh in 1887, who determined the so-called Rayleigh waves. This work was later extended 
for propagation of waves in elastic plates by Lord Rayleigh (1889) himself and Lamb (1889). 
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Here we extend the Rayleigh-Lamb problem and account for incremental motions superimposed 
on finite deformations in dielectric media, and investigate how these are influenced by the presence 
of external electric field and pre- stretch. 

When subjected to an electric field electroactive polymers (EAPs) deform and both their me- 
chanical and electrical properties are modified. In contrast to piezoelectric ceramics, DEs are capa- 
ble of undergoing large deformations, a property that entitled them the name "artificial muscles". 
Moreover, while in piezoelectricity the electromechanical coupling is linear, in DEs the mechanical 
field depends quadratically on the electric fields. A proper theory, which accounts for the afore- 
mentioned coupling and captures the abihty of the material to undergo finite strains is therefore 
required. The foundations of this nonlinear electroelastic theory are summarized in the pioneering 
works of Toupin (1956) and Eringen (1963) for the static case. These contributions were later ex- 
tended by Toupin (1963) to account for the dynamics of these elastic dielectrics. A comprehensive 
summary can be found in monographs by Eringen and Maugin (1990) and Kovetz (2000). Due to 
the development of new materials that admit this coupled behavior, thus branching toward a win- 
dow of new applications (e.g., Pelrine et al., 2000, Bar-Cohen, 2002 and Carpi and Smela, 2009), 
the interest in these electroelastodynamic theories revived recently and the coupled electromechan- 
ical theory was revisited recently (e.g., McMeeking and Landis, 2005; Dorfmann and Ogden, 2005; 
Ericksen, 2007; Suo et al., 2008). The foregoing works and their extension to the domain of soft 
dielectric composites by deBotton et al. (2007) and Bertoldi and Gei (2011), differ in their consti- 
tutive formulations, the choice of the independent variables, the resultant electrostatic stress-like 
tensors and electric body-like forces. For a review of the diverse formulations of the constitutive 
laws and governing equations the reader is referred to Bustamante et al. (2009). Among the various 
approaches we recall the formulation proposed by Dorfmann and Ogden (2005) for the static case, 
and its extension to dynamics by Dorfmann and Ogden (2010). In these works the concept of 'total' 
stress tensor that is derived from a 'total' or 'augmented' energy-density function was employed. 
In this work we foUow the framework proposed by the latter. 

In contrast with the relatively large body of theoretical works that are available, only a few 
boundary-value problems (BVP) were solved in the context of the dynamic behavior of electrome- 
chanically coupled EAPs. One of the first contributions in the field of piezoelectricity was made by 
Tiersten (1963a), who examined the thickness vibrations of an infinite piezoelectric plate induced 
by alternating voltage at the surface electrodes, and later solved the corresponding problem of wave 
propagation (Tiersten, 1963b). To the best of the authors knowledge, while for piezoelectric media 
solutions for BVP accounting for the influence of pre-strain and bias field are available (see review 
article by Yang and Hu, 2004), analogue developments in the context of dielectrics were considered 
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by Mockensturm and Goulbourne (2006) and Zhu et al. (2010), who examined the dynamic behav- 
ior of dielectric elastomer balloons, and by Dorfmann and Ogden (2010) who studied the problem 
of propagation of Rayleigh surface waves in a dielectric half-space. Herein we continue along the 
path of the latter contribution and consider the extension of the Rayleigh-Lamb wave propagation 
problem to finitely deformed dielectric layers subjected to coupled electromechanical loading. 

The work is composed as follows. In section 2 the theory of finite electroelastodynamics is sum- 
marized. The equations for incremental motions superimposed on finite deformations are outlined 
in section 3. Specific finite deformations that correspond to three different loading configurations 
are considered in section 4 for a layer whose behavior is characterized by a particular augmented 
energy-density function (AEDF), namely the incompressible dielectric neo-Hookean model (DH). 
Next, the extension of the Rayleigh-Lamb dispersion relation for a DH layer is introduced in sec- 
tion 5. An analysis of the dispersion relation is carried out in section 6, and the effects of the bias 
electrostatic field and pre-stretch are investigated for the three loading configurations. The main 
conclusions and observations are summarized in section 7. 



2. Finite electroelasticity 

Let ;t; : X ^ C M"^ describe the motion of a material point X from a reference configu- 
ration of a body Q.o, with a boundary dQ.o, to its current configuration Q, with boundary dQ., by 
x = X (X, t), where is a time interval. The domain of the space surrounding the body is ]R^\f2 and 
is assumed to be vacuum. The corresponding velocity and acceleration are denoted hy \ = x t ^i^d 
3 = Z ;n respectively, while the deformation gradient is F = ^ = Vx;^, and where due to the ma- 
terial impenetrability J = det (F) > 0. Vectors between two infinitesimally close points are related 
through dx = FdX, whereas area elements are transformed via Nanson's formula NdA = jF^nda. 
The volume ratio between an infinitesimal volume element dv in the deformed configuration, and 
its counterpart in the reference dV is given by dv = JdV. As measures of the deformation the right 
and left Cauchy-Green strain tensors C = F^F and b = FF^ are used. 

Let e denote the electric field in the current configuration. Commonly the electric field is given 
by means of a gradient of a scalar field, namely the electrostatic potential. The induced electric 
displacement field d is related to the electric field in free space via the vacuum permittivity Cq such 
that d = Coe. In dielectric media an appropriate constitutive law specifies the relationship between 
these fields. Generally, this connection can be non-linear and anisotropic. 

The balance of linear momentum is 

V.c7 = pa, (2.1) 
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where o is the 'total' stress tensor, and p is the material mass density. The balance of angular 
momentum implies that a is symmetric. Note that a consists of both mechanical and electrical 
contributions, such that the traction t on a deformed area element can be written as en where n is 
unit vector normal to dQ.. On the boundary of the material we postulate a separation of the traction 
into the sum of a mechanical traction t,„ which is a prescribed data, and an electrical traction te 
which is induced by the external electric field. 

Assuming no free body charge (ideal dielectric). Gauss' law reads 

V-d = 0. (2.2) 

Under a quasi-electrostatic approximation, appropriate when for the same frequency the length of 
the waves under consideration are shorter than the electromagnetic waves, Faraday's law states that 
the electric field is curl-free, i.e., 

V X e = 0, (2.3) 

thus enabling the usage of the aforementioned electrostatic potential. 

Taking into account fields outside the material, which henceforth will be identified by a star 
superscript, the following jump conditions should be satisfied across dQ., namely 



[[a]] n = t„ 



[[d]]-n=-w„ [[e]]xn = 0, 



(2.4) 



where Wg is the surface charge density, and the notation [[•]] = (•) — (•)* is used for the difference 
between fields inside and outside of the material. The outer fields are related by 



Coe" 



£0 



-* "*)I 



e -e" 



(2.5) 
(2.6) 



where I is the identity tensor. Herein we identify the electrical traction tg as the consequence of 
the external stress a*, namely the Maxwell stress, such that tg = o*n. In the surrounding space 
outside the material d* and e* must satisfy Eqs. (2.2)-(2.3), which reduce to Laplace equation of 
the electrostatic potential. As a consequence the Maxwell stress is divergence-free. 

The foregoing balance and jump equations can be recast in a Lagrangian formulation with the 
appropriate pull-back operations. Specifically, we have that 



P = 7aF 



-T 



D = 7F Id 



E = F^e, 



(2.7) 
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for the 'total' first Piola-Kirchhoff stress, Lagrangian electric displacement and electric field, re- 
spectively (e.g., Dorfmann and Ogden, 2005). The corresponding balance equations are 

Vx-P = PLa, Vx-D = 0, VxxE = 0, (2.8) 

where = Jp is the density of the material in the reference configuration. The jump conditions 
across the boundary dQ.o read 

[[P]]N = tM, [[D]]-N = -W£, [[E]]xN = 0, (2.9) 

where tMdA = t^da, w^dA = Wgda and N is a unit outward normal to DUq. 

Following Dorfmann and Ogden (2005), the 'total' first Piola-Kirchhoff stress and the Lagrangian 
electric field are given in terms of an augmented energy-density function ^ (AEDF) with the inde- 
pendent variables F and D, such that 

For an incompressible material a Lagrange multiplier p is introduced, which is a workless reaction 
to the kinematic constraint such that 

P=--;.F-T. (2.11) 
The latter can be determined only from the equilibrium equations and the boundary conditions. 

3. Small fields superimposed on finite deformations 

Recent experiments reveal how pre- stretching of dielectric elastomers enhances properties such 
as the actuation strain (Pelrine et al., 2000) and breakdown strength (Plante and Dubowsky, 2006). 
A thorough investigation of the subject was recently done by Kofod (2008). Motivated by these 
findings, we address the response of pre-stretched dielectric elastomers to incremental deforma- 
tions. 

Following the formulation given in Dorfmann and Ogden (2010), we consider infinitesimal time- 
dependent elastic displacement and electric displacement increments x = (X,?) and D (X,?), re- 
spectively, superimposed on the foregoing static configuration reached via x (X) . Herein and 
throughout this work a superposed dot will denote the incremental quantities. The corresponding 
balance laws can be formulated in terms of Eulerian quantities, namely 

V-E = px„, V-d = 0, Vxe = 0, (3.1) 

such that E d and e are, respectively, the push-forwards of the increments in the total first Piola- 
Kirchhoff stress, the Lagrangian electric displacement and the electric field determined via the 
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inverse of the transformations given in Eq. (2.7). Next, upon linearization, the incremental consti- 
tutive equations for an incompressible material are 

Z = 'tg'h + ph^ -pl + M, (3.2) 
e = ^^h + £^d, (3.3) 

where (e^^h)^ = I^ij^hij. For a compressible material Eq. (3.2) is to be taken with p = p = 0. 
Herein h = Vx is the displacement gradient and 

■^ij = -^Pai^^pl-^OaP^ ^ijk= FjaFpi^^Oiaji, %jkl = jFjaFi^'^Qiak^^ (3-4) 

are the push-forwards of the referential electric, electroelastic, and elasticity tensors, respectively. 
The latter are defined by 

^2ij, ^2vj/ ^2vj/ 

Similarly, the incremental outer fields are 

d* = Coe*, (3.6) 
a* = £0 [e*(g)e* + e*®e*- (e*-e*)l] , (3.7) 

where d* and e* are to satisfy Eqs. (2.2)-(2.3), hence C7* is identically divergence-free. 

Utilizing several kinematic relations (see Dorfmann and Ogden, 2010) the corresponding jump 
conditions are 

[E-c7* + c7*h'^-(V-x)(7*]n = tm, (3.8) 
[d-d*-(V-x)d'' + hd*] -n = -w,, (3.9) 
[e-e*-hV] xn = 0, (3.10) 

where the quantities tm and Wg are defined by t„,da = tMdA and vv^da = vv^dA. 

4. Incompressible dielectric neo-Hookean layer subjected to transverse bias 

FIELD 

We recall that the constitutive laws in terms of the AEDF has the deformation gradient and 
the electric displacement field as its independent variables. Whenever the material is isotropic 
its symmetry can be exploited, and the dependency of the AEDF on F and D is reduced to the 
appropriate invariants of C and D ® D, namely 

/i=tr(C)=C:I, /2 = i(/i2-C:C), I3 = dct{C) = J^, (4.1) 
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and 

/4, = tr(D®D), /5, = C:(D®D), /g, = : (D®D) . (4.2) 
Thus, we consider the AEDF 

^D// {hJse) = I (/i - 3) + ^/5„ (4.3) 

where /i is the shear modulus and £ is the material permittivity that equals the vacuum permittivity 
Co times the material relative permittivity e^. This model recovers a neo-Hookean behavior in 
a purely elastic case, as well as the simple isotropic linear relation d = ee between the current 
displacement and electric field. Henceforth we refer to Eq. (4.3) as the incompressible dielectric 
neo-Hookean model (DH). Throughout this work the material behavior is assumed to be governed 
by the DH model. The corresponding total stress is 

a = ^b + ^d(g)d-pl, (4.4) 

where, in component form, the constitutive tensors , and are 

= ^5//, ^ijk = I {dikdj + didjk) , %jk! = IJ-Sikbji + ^dikdjdi. (4.5) 

Consider a Cartesian coordinate system with unit vectors ii,i2 and i3 along the Xi,X2 and 
axes, respectively. Let a layer governed by the DH model be infinitely long along the xi-axis with 
a thickness 2h along the X2-axis. Assuming a plane-strain configuration there is no deformation 
along the X3-axis and the fields being independent of X3. An electric displacement field is applied 
to the layer along its thickness. In the current configuration this field is given by d = ^2^2^ and 
the corresponding electric field equals the voltage difference between the electrodes divided by 
2h. The resultant deformation in terms of a homogeneous diagonal deformation gradient F = 
diag [A, \ 1] , is related to the stress via 

aii=juA2-p, 022 = ji + ^d^- p. (4.6) 
Specialization of the in-plane components of Eq. (3.5) for the considered deformation gives 

Ml = =2^22 = p (4.7) 

=^121 = =5^211 = ^=^222 = ^^;?2, (4.8) 
^1111 = ^2121 =M', (4.9) 
^1212 = "^2222 = + -d2- (4.10) 
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Figure 4.1. Illustration of (a) path A: expansion-free loading, (b) path B: initially 
pre-stretched layer, and (c) path C: a layer immersed in electric field. 

Wave propagation and vibrations of the layer are addressed for the following three fundamental 
loading paths. 

Path A: expansion-free finite loading. The layer having its top and bottom surfaces = ±/? 
coated with soft electrodes and mechanically traction free (t^ = 0). It is also free to expand along 
the xi-axis (Fig. 4. la). From the traction-free boundary conditions simple connections between the 
stretch, the pressure and the Lagrangian electric displacement field D = Ad are obtained, namely 

X = {l+D^f\ p = ^{l+D^f\ (4.11) 

where D = D2/y/ll£. Note that due to the symmetry of the problem the electric fields along with 
the Maxwell stress outside the layer vanish. The surface charge we is evaluated via Eq. (2.9), 
recalhng that d* = 0. 

Path B: initially pre-stretched layer. Once again the surfaces X2 = ±h are coated with electrodes 
and are free of mechanical traction. The layer is first pre-stretched along the xi-axis to A = A, 
and then clamped (Fig. 4.1b). Subsequently, an electric displacement field is applied, while the 
longitudinal stretch is kept constant. The pressure is given in terms of the pre-stretch and the 
electric displacement field as 

p = ^Q^ + cpy (4.12) 

where d = ^J\X£-- The resultant pre-stress is obtained via Eq. (4.6). Again, the electric fields 
and the Maxwell stress outside the layer vanish. The surface charge We is evaluated via Eq. (2.4), 
having that d* = 0. 

Path C: a layer immersed in electric field. The layer is immersed in a pre-existing electric field. 
The surfaces are free of mechanical traction. The absence of the electrodes on the surfaces im- 
plies that We = 0, where the pressure and the stretch are determined via the corresponding jump 



conditions to yield 



X = {\-b\er-l))"\ p = ^{l+X'). (4.13) 



We note that when deriving the solution it was assumed the electric field is homogeneous outside 
the layer, neglecting fringing effects. Thus, this solution is by no means exact, and serves only as a 
qualitative result. Eq. (4.13a) reveals that the layer contracts along the xi-axis and expands along 
the direction of the electric displacement, contrary to the situation in path A. 

5. Wave propagation in dielectric layers 

We determine next the mechanical and electrical waves due to harmonic excitation along the 
xi-axis superimposed on the deformed configurations described previously. Following Dorfmann 
and Ogden (2010), the incremental equations of motion (3.1a) read 

[jiX^ +p)hn,i -p,i+ (|l"^i^2^ hi2,2+ph2i,2 + ^d2di,2 = pxi,tt, (5.1) 

{llX^ + p) h2lA + Phl2,\ + ^d2di^i + + + p^ h222 - P,2 + '2^d2d2,2 = pi2,K, (5.2) 

while the incremental Faraday's equation (3.1c) is 

^d2ihi2,2 + h2i,2) + ^di2-2^d2h22,i-^d2,i =0. (5.3) 

The motion has to satisfy the incompressibility constraint and Gauss' equation, namely 

hn+h22 = 0, ^1,1+^2,2 = 0, (5.4) 

which motivate the use of stream functions ^ {xi,X2j) and cp {xi,X2j) such that 

^1=^2, i2 = -fi, di = (P,2, d2 = -(p,i. (5.5) 

Normalization by /iA^, followed by differentiation of Eqs. (5.1) and (5.2) with respect to X2 andxi, 
respectively, and then subtraction of the latter form the former results in the coupled equation 
\2\ 1 , n2 



/ 1+Z)2\ 1+^^, D . , 1 p 

0,1111+1^1 + —^ j (^,1122+ ^4 0,2222 + ^3^^ (<?>,112 + 9,222) = p - (0,1 1 + 0,22) ^„ 



(5.6) 



After multiplying Eq. (5.3) by A ^ / ^ it may be rewritten in the form 

V r*- 



^(0,112 + 0,222) + ^(9,11 + 9,22) = 0. (5.7) 

The last two equations constitute a system of two coupled equations for the stream functions 
and (p. A solution is feasible once the current stretch A and the dimensionless nominal electric 



displacement D are prescribed in addition to the density p and shear modulus /i of the layer. We 
recall that while for the paths A and C the quantities D and A are uniquely related via Eqs. (4.1 la) 
and (4.13a), respectively, they are independent for path B. Also note that the quantity 1 / y/JIs can 
be integrated into the amplitude of the stream function cp. 

Next we assume a solution with time-dependency in the form e^'^' and a periodicity along the 
xi-axis in the form e'^^^, where O) is the angular frequency and k is the associated wavenumber, 
such that the wave velocity is c = CO /k. Thus, travehng waves solutions in the form 

are sought, where the appropriate values of q are to be found from the governing equations. In- 
sertion of Eq. (5.8) into Eqs. (5.6)-(5.7) yields the following system of two linear homogeneous 
equations in A and B, namely 

{q'-l)(^^c' + q'^-^-X'y+{q'-l)qlDB = 0, (5.9) 

q{q^-\)DA+{q^-l)?iB = 0, (5.10) 

where B = B/y/Jie. A non-trivial solution is possible for values of q for which the determinant 
of the coefficients of A and B vanish. This leads to a bi-cubic equation in q. The six roots in- 
clude two repeated ones, namely qi = qj = I and ^3 = ^4 = — 1, along with roots qs = —qe — 

An interesting case arises if = A^^, that is the square of the speed of a bulk shear wave in an 
isotropic elastic medium subjected to a finite deformation propagating along the axis with principal 
stretch A . In this case qs = qe = 0, and the corresponding solution is associated with propagation 
of waves parallel to the boundary along the xi-axis. As observed by Dorfmann and Ogden (2010) 
for the half-space case, the associated stream functions vanish as expected since an in-plane bulk 
shear wave cannot propagate in this manner. 

Usually, the general solutions for ^ and (p are obtained as a linear combination of the roots in 
the form of 

6 6 
(j) = (p = £ (5.11) 

n=l n=l 

However, a somewhat different structure of the solution is required for the DH model. For the 
general case when the constitutive law involves l4e, a full coupling between the mechanical and 
electric equations is expected (see Dorfmann and Ogden, 2010), and hence dependency between 
An and 5„ exists for all n. However, for the specific constitutive behavior considered here we 
find that for the roots {qn}, n= 1, ..,4, the coefficients associated with Ai, ...,A4, 5i, ...,54 in the 
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counter diagonal of the coefficients matrix vanishes. Mathematically, this implies that the constants 
Ai, ...,A4 and 5i, ...,84 are independent. Physically, this hints at the possibility of propagation of 
elastic waves without excitation of the accompanying electric waves. Thus, the associated solutions 
for and cp of these roots are not to be taken as multiple. The coupling between the coefficients 
still holds when ^5 and q(, are considered, giving a relation between A5 and ^5 , and between Ag and 
B^. This relation is obtained by substituting ^5 and ^5 into Eq. (5.9), or equivalently into Eq. (5.10). 
Hence, the solutions sought for the DH layer are 

(2=1,3,5,6 77=1, 3, 5, 6 

where the six unknowns are Ai,A3,A5,A6,5i and ^3, while and B(, are functions of A5 and Ag. 

We note that once the layer is excited the symmetry in paths A and B is broken. Accordingly, 
for all paths the incremental exterior fields must be accounted for. Motivated by the need to satisfy 
a decaying condition at X2 — > ±0° along with Laplace's equation outside the material, we consider 
the following stream functions 

V/* = Ci/yte"^^2e*(-'^'"''') atJC2 >/z, (5.13) 
t}* = C2/fe^''2g^/:(^i-rt) atx2 < -/z, (5.14) 

such that the electric field components are given by 

{-xir*, X2>h \ X2>h 

' ' ^2= ' • (5-15) 

-1^*1 X2 < -h {-^*2 ^2 < -h 

In summary, we end up with a set of eight constants to be determined, namely {A„} , n = 1 , 3, 5, 6, 
{Bn},n =1,3 and Ci, C2 via the appropriate jump conditions. Henceforth we distinguish the 
solutions of paths A and B from the one along path C. 



Paths A and B. The Maxwell stress, and as a consequence of Eq. (3.7), its increment vanish 
outside the layer (a* = C7* = 0). Further, the surfaces remain traction-free and hence the jump in 
the stress according to Eq. (3.8) becomes 

L22 = 0atX2 = ±h, (5.16) 

Z12 = 0atX2 = ±h, (5.17) 
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where 

^22 = - {j2 + l'^2+P^ ^^n~p-2^d2(p,u (5.18) 

2:12= (^j2 + ^4^^^22-p^,n + ^d2(p,2. (5.19) 
The solution sought for the incremental pressure is 

p{xi,X2,t) = k (^1/^1^2+^3/^3x2^ g*(xi-cO^ (5 20) 

where Pi and ^3 are determined via Eqs. (5.1)-(5.2). 

As the electrodes surface charge is fixed we have that vvg = 0, and the corresponding jump 
condition in Eq. (3.9) becomes 

d2 — d2 = 0atX2 = ±/?, (5.21) 

where d^ = £0^2- The jump in Eq. (3.10) is 

e\-ei=0atX2 = ±h, (5.22) 

where ei = ^^2 (0,22 — 0,ii) + ^<P,2. Note that while the increment in the Maxwell stress is zero as 
discussed before, Eqs. (5.21)-(5.22) state that d* and e* do not vanish. 



Path C. Along this path the Maxwell stress and its increment do not vanish outside the layer, 
yielding the following form of Eq. (3.8) 

E22-C7*22 + Cy2*2^22 = at X2 = ±h, (5.23) 

J^i2-<7*i2 + <7n/i2i = 0atX2 = ±h, (5.24) 

where 

= -c^22 = -;^. c^*22 = £0^2^' a*i2 = £oete^. (5.25) 

2£o 

The expressions for E22 and E12 are given in Eqs. (5.18)-(5.19). The external electric fields do not 
vanish in this path, and hence Eq. (3.9) assumes the form 

dl-d2 + d2h22 = 0atX2 = ±/j, (5.26) 

where = '^2- The remaining jump in Eq. (3.10) is 

e\-ei+e^h2i=0atX2 = ±h. (5.27) 
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Eqs. (5.16)-(5.22) andEqs. (5.23)-(5.27) complete the necessary set of eight boundary conditions 
for the eight unknowns for paths A and B, and path C, respectively. Furthermore, they constitute 
a set of linear homogeneous equations in the unknowns, hence non-trivial solutions exist when the 
determinant of the coefficients matrix vanishes. This, in turn, is the extension of the well-known 
Rayleigh-Lamb transcendental equation for a purely elastic layer, and Tiersten's result for a piezo- 
electric plate (Tiersten, 1963a) within the framework of infinitesimal deformations. Of practical 
interest is the dependency of the dispersion relation, usually given in terms of the frequency spec- 
trum, on the bias field and the finite deformation. These are discussed in the next section 

6. Analysis of the dispersion relation 

The working scheme of this section can be formulated as follows: Given an angular frequency 
of excitation CO, we determine the associated wavenumbers {k^} satisfying the generalized tran- 
scendental equation, determine the velocities of the propagating waves Cm = (O/ km, and examine 
how these vary as functions of the bias field and pre-stretch. Specifically, for paths A and C the 
dependency on the bias field is examined for a few representative values of D. For path B, having 
the deformed configuration fixed we find it advantages to explore the dependency in terms of d. 
Recalling that along this path the pre-stretch and the electric displacement field are independent, 
the dependency of the dispersion relation on a few representative values of the pre-stretch X is also 
examined. Of special interest is the identification and evaluation of the fundamental modes, i.e., 
those modes having finite velocity in the limit of long waves (e.g., Lutianov and Rogerson, 2010). 

We find it useful to represent the results in terms of the dimensionless quantities k = hk, c = c/cb 
and (b = ck, where cb = \/l^/p is the bulk shear wave velocity in an isotropic elastic material. 
Further, our numerical investigation implies that the value of the parameter has only minor a 
influence on the frequency spectrum due to the choices D = D2/ yJJiE and J = J2/ yJV^ which are 
normalized by the material permittivity. Accordingly, as a representative value we chose = 3. 

In their fundamental works Lord Rayleigh (1889) and Lamb (1889) had shown that in the purely 
elastic case a decomposition of the mechanical vibrations can be made into symmetric (exten- 
sional) and antisymmetric (flexural) modes with respect to the mid-plane of the layer. Interestingly, 
a similar decomposition of the mechanical displacement is possible for the coupled problem. Inter- 
estingly, it turns out that the symmetry of the electrical displacements is reversed in the sense that 
whenever the mechanical displacements are symmetric, the electrical displacements are antisym- 
metric, and vice versa. To distinguish between the two types of modes we employ the following 
procedure. For the antisymmetric modes we subtract Eq. (5.23) evaluated at X2 = —h from its value 
at X2 = h, and similarly with Eq. (5.27). Subsequently, we add Eq. (5.24) evaluated at X2 = —h to 
its value at X2 = h, and similarly with Eq. (5.26). We end up with a system of four homogeneous 
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equations in terms of the four constants 

A+=Ai+A3, A+=A5+A6, B =B3~Bi, C+=Ci+C2. (6.1) 

The dispersion relation derived from the determinant of this system provides the speed of propa- 
gation of waves associated with the antisymmetric modes. The symmetric modes can be derived 
by changing the foregoing operations. Thus, to add Eq. (5.23) evaluated at X2 = —h to its value at 
X2 = h, and similarly with Eq. (5.27). Analogously, to subtract Eq. (5.24) evaluated atX2 = —h from 
its value at X2 = h, and also with Eq. (5.26). These operations lead to system of four homogeneous 
equations in terms of the four constants 

A-=Ai-A3, A^=A^-Ae, B+=B3+Bi, C =Ci-C2. (6.2) 

Path A: expansion-free finite loading. Fig. 6.1a displays the normalized velocity c as a function 
of the normalized wavenumber k for a few values of D. Herein and henceforth the continuous and 
dashed curves correspond to the symmetric and antisymmetric modes, respectively. The markless 
curve, and the curves with triangle and circle marks correspond to D = 0, 1 and 2, respectively. 
These values of the biased fields induce the principal stretches A = 1,1.19 and 1.49, respectively. 

The case D = corresponds to the purely elastic problem. In this case the generalized solution 
recovers the classic Rayleigh-Lamb dispersion relation, as it should. Note that one can perceive 
variation in k as variation in the plate thickness for a fixed wavenumber. Thus, as k increases 
the plate geometry becomes similar to that of a half-space. Accordingly, the symmetric and the 
antisymmetric branches become closer, such that in the limit of ^ — )■ oo they coincide and attain 
the surface wave velocity. In particular, the fundamental modes propagate at a velocity c = 0.955 
in the limit of short waves (large k), which is the propagation velocity of Rayleigh surface waves. 
Further, the fundamental symmetric mode attains the well known result c = 2 in the limit of long 
waves, or, equivalently, thin layers (small k). 

As D increases a monotonous rise of the velocity in the limit of long waves for the symmetric 
mode is revealed. The velocity at this limit can be derived analytically upon taking the first term 
in the Taylor series expansion of the symmetric dispersion relation in the neighborhood of ^ = 0. 
Equating it to zero yields an explicit expression for the velocity in the limit of long waves, namely 

c = 2{l+D^y^^ = 2X. (6.3) 

A monotonous rise of the velocity in the limit of short waves is observed too. We stress that 
the asymptotic values of the dimensionless speed are in agreement with the corresponding results 
obtained with an analysis of surface waves. 
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Figure 6.1. Path A: the normalized velocity c (a) and angular frequency d) (b) as 
functions of the normalized wavenumber k for the expansion-free finite loading path. 
The continuous and dashed curves correspond to the symmetric and antisymmetric 
parts, respectively. The markless curve, and the curves with triangle and circle 
marks correspond to D = 0, 1 and 2, respectively. 

We find that all antisymmetric branches emerge at the origin, so that no loss of stability in terms 
of diffuse mode is identified, i.e., a vanishing c for non-zero wavenumbers. For this loading path a 
possible instability could be macroscopic instability associated with the loss of positive definiteness 
of the tangent constitutive operator (Zhao and Suo, 2007, Bertoldi and Gei, 2011) that however 
cannot occur for a DH layer with vanishing longitudinal stress. 

Fig. 6.1b displays the normalized angular frequency G) as function of the normalized wavenum- 
ber k. We point out that the slope of the curves remains positive. Consequently, a unique wavenum- 
ber k is associated with each value of CO along a specific curve. In other words, exciting the layer 
along this path with a certain frequency will give rise to a single wavelength of each fundamental 
mode, and increasing values of OJ will result in shorter waves. 

Path B: Pre-stretched layer. Three representative types of pre-stretch where chosen, namely a 
tensile pre-stretch with A = 1.5, no pre-stretch A = 1 (the layer is clamped in its the reference 
configuration), and a pre-compression with A = 0.8. We recall that here we find it useful to use 
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Figure 6.2. Path B: the normalized velocity c (a) and angular frequency &> (b) 
as functions of the normalized wavenumber k for for A=1.5. The continuous and 
dashed curves correspond to the symmetric and antisymmetric parts, respectively. 
The markless curve, and the curves with triangle, circle and square marks corre- 
spond to 6? = 0, 1, 2 and 2.24, respectively. 



the current quantity d = y/JiE as a measure for the bias field, as A is held fixed. Figs. 6.2a, 
6.3a and 6.4a display the normalized velocity c as a function of the normalized wavenumber k for 
A = 1.5, 1 and 0.8, respectively. The continuous and dashed curves correspond to the symmetric 
and antisymmetric parts, respectively. The markless curve, and the curves with triangle and circle 
marks correspond to J = 0, 1 and 2, respectively. The curves with square marks correspond to 
d = 2.24, 1.75 and 1.44 in Figs. 6.2-6 A, respectively. 

With regard to the symmetric branches, the velocity in the limit of long waves increases with 
values of d, in a manner similar to the one observed for path A. Here again we take the first term in 
the Taylor series expansion of the symmetric dispersion relation in the neighborhood of ^ = and 
equate it to zero. The resulting explicit expression for the velocity in this limit is 

c = ^3 (j2 + A-2) +A2. (6.4) 
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Figure 6.3. Path B: the normalized velocity c (a) and angular frequency (b (b) as 
functions of the normalized wavenumber k for A=l. The continuous and dashed 
curves correspond to the symmetric and antisymmetric parts, respectively. The 
markless curve, and the curves with triangle, square and circle marks correspond 
tO6? = 0,l,1.75 and 2, respectively. 

Interestingly, a reversed trend is revealed in the limit of short waves, such that the surface wave 
velocity decreases monotonically with d, in contrast to the situation in path A. The difference stems 
from the emergence of compressive stress in path B, as the layer cannot elongate with increasing 
values of J due to the fixed boundary. The surface wave velocity decays until it reaches a vanishing 
value associated with a loss of surface stability. This is attained when a threshold value dth is 
applied, such that 

lim c = 0. (6.5) 

Thus dfh can be perceived as the value of J at which a half-space losses its stability. This threshold 
value depends on the pre-stretch such that it increases monotonically with A, revealing a stabilizing 
effect of the pre-stretch. Specifically, for A = 1.5, 1 and 0.8 we calculate numerically the values 
dth = 2.245, 1.753 and 1.453, respectively. When d > dth the velocity vanishes at a finite cutoff 
wavenumber which we denote by ^co™- In a way Ic/o'" defines a critical wavelength, beyond which 
no propagation of waves is exhibited. Further increase of d will result in smaller values of ^co'"- For 
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Figure 6.4. Path B: the normalized velocity c (a) and angular frequency (b (b) as 
functions of the normalized wavenumber k for A=0.8. The continuous and dashed 
curves correspond to the symmetric and antisymmetric parts, respectively. The 
markless curve, and the curves with triangle, square and circle marks correspond 
to 6? = 0, 1, 1.44 and 2, respectively. 

completeness, we note that our numerical investigation shows that for the case d = the solution 
recovers Blot's result for loss of stability under compression of a neo-Hookean half-space at a 
pre-compression of A = 0.544. 

Turning to the antisymmetric branches, we observe an essentially different evolution from the 
one in path A. Thus, in this case the branches emerge from various points. We denote the wavenum- 
ber at which the branch emerges by kco'\ In a way k^"'^ defines a critical wavelength beneath which 
there is no stable propagation of longer waves. The value of A;™ increases when d is increased, up 
to the threshold value dtu, at which tends to infinity. We stress that dtu is the one satisfying 
Eq. (6.5) as expected, since in the limit of half-space the two branches coincide. 

Figs. 6.2b, 6.3b, and 6.4b display the dispersion relation in terms of the normalized frequency G) 
as a function of the normalized wavenumber k, with identical legend to the one previously used. It is 
observed how, unlike the corresponding curves in path A, for some values of J the symmetric curves 
are not monotonous. In some cases the slope of the curves is changing from positive to negative 
and to positive once again. This corresponds to a situation where for the same frequency there 
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Figure 6.5. Path C: the normalized velocity c (a) and angular frequency (b (b) as 
functions of the normalized wavenumber k for a layer immersed in an electric field. 
The continuous and dashed curves correspond to the symmetric and antisymmetric 
parts, respectively. The markless curve, and the curves with triangle, circle and 
square marks correspond io D = 0,0.5,0.632 and 0.65, respectively. 

are more than one wavenumber satisfying the dispersion relation for the fundamental symmetric 
mode, and subsequently the pertained different wavelengths propagate with different velocities. If 
the frequency is to be taken as the independent variable, this suggests a peculiar phenomenon in 
which as the frequency is increased there is a possibility for the appearance of shorter and slower 
waves that propagate along with the "primary" waves. 

Path C: A layer immersed in an electric fields. We recall that Eq. (4.13a) states that the layer 
expands along the direction of the electric displacement, in a manner opposite to the one revealed 
in paths A and B. Further, the stretch depends strongly on the material permittivity, even when the 
normalized quantity D is used. For the choice = 3 the maximal value of admissible D is 1/ \/2, 
as X vanishes for this value in view of Eq. (4.13a). Consequently, the values of D considered in 
path A cannot be used, and even values of D that result in the same values of X are not feasible as 
here A < 1, whereas along path A A > 1. Nonetheless, attempting to investigate cases analogue 
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to those addressed in path A, we choose values of D that yield values of max {A, 1/A} similar to 
those considered in path A. 

Fig. 6.5a displays the normalized velocity c as a function of the normalized wavenumber k for 
D = 0,0.5,0.632 and 0.65. These values correspond to principle maximal stretche ratios 1, 1.2, 1.5 
and 1.59, respectively. The continuous and dashed curves correspond to the symmetric and an- 
tisymmetric solutions, respectively. The markless curve, and the curves with triangle, circle and 
square marks correspond to D = 0,0.5,0.632 and 0.65, respectively. Once again, the symmetric 
modes exhibit a rise in the velocity in the limit of long waves with increasing values of D. In a man- 
ner similar to path B, the trend is reversed in the short waves limit, and the surface wave velocity 
decreases monotonically as D increases. This is a consequence of the resultant compressive stress 
that can be determined by substitution of Eq. (5.25a) into Eq. (2.4a) for the jump in the stress. Thus, 
the surface wave velocity decays as a function of D, until it reaches a vanishing value associated 
with a loss of surface stability, i.e., satisfying Eq. (6.5). Once again we denote the threshold value 
of D for which the latter holds by D^/,. In this numerical example we find that Df^ = 0.636. As in 
Path B, when D > Dth the velocity of the symmetric modes vanishes at a finite cutoff wavenumber 
ka"^ beyond which no shorter waves are exhibited. 

The behavior of the antisymmetric modes is also reminiscent of the one observed for path B. 
Specifically, it is observed how the curves emerge at different cutoff wavenumbers k""^' defin- 
ing critical wavelength beneath which there is no stable propagation of longer waves. The cutoff 
wavenumber increases monotonically with values of D, and reaches infinity when D attains the 
aforementioned value Dth- Subsequently, when D > Df^ there is no stable antisymmetric wave 
propagation. 

Fig. 6.5b displays the dispersion relation in terms of the normalized frequency O) as function 
of the normalized wavenumber k. Trends similar to the ones observed in path B are revealed. In 
particular, the monotonicity is lost when D > 0.584, where the positive slop of the curve becomes 
negative and then positive again. Here again, this implies that imposing certain frequencies may 
yield more than one wavelength of the fundamental symmetric mode satisfying the dispersion rela- 
tion. In turn, this suggests that when the frequency is increased there is a possibility of simultaneous 
different phase velocities. 

7. Concluding remarks 

In view of the emergence of modem EAPs capable of large deformations in response to elec- 
tric excitation we examine the topic of small waves propagation superposed on finitely deformed 
dielectric elastomer layers. Specifically, we considered a layer whose electromechanical behavior 
is characterized by the DH model when subjected to a coaxial finite deformation. Three different 
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loading paths which result in the aforementioned finite configurations were addressed. In path A the 
layer is free to expand along its longitudinal direction and deformed due to the electric excitation 
along its thickness. In path B the layer is initially pre- stretched, then clamped, and subsequently 
the electric field is employed. In path C the layer is immersed in a pre-existing electric field. 

Following Dorfmann and Ogden (2010) small perturbations on top of the homogeneous finite 
deformation are considered. The coupled equations of incremental electrodynamics along with 
the pertained boundary conditions yielded the desired dispersion relation. For the DH model this 
required a non-conventional presentation of the solution of the governing equations. It was shown 
that a separation of the mechanical and electric waves to symmetric and antisymmetric modes with 
respect to the mid-plane of the layer is feasible. We find that the symmetric mechanical waves are 
accompanied by antisymmetric electric waves and vice versa. 

Numerical investigation of the fundamental modes was conducted to examine the influence of 
the bias field along the three loading paths. For path B we also examined the effect of pre-stretch. 
In path A the velocity of both short waves (half-space, large k) and long waves (thin plate, small k) 
is increased when the electric displacement field is enhanced. While in the limit of long waves a 
similar rise in the velocity is observed in path B, the velocity decreases as function of the electric 
displacement field in the limit of short waves. The existence of a threshold values of J at which the 
layer losses stability is observed. The pre-stretch effect is shown to be stabilizing, as higher values 
of d are needed to reach the onset of instability. Dispersion relations similar to the ones obtained 
in path B were observed in path C. We finally stress that along all paths there is a marked influence 
of the bias electric field and pre-stretch on the propagation speed. This phenomenon lends itself as 
a possible control mechanism for the speed of the waves, or even for filtering specific wavelengths 
with suitable adjustment of the bias field and pre-stretch. 
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